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ABSTRACT 

NASA’s Kepler mission has provided several thousand transiting planet candidates during the four 
years of its nominal mission, yet only a small subset of these candidates have been confirmed as true 
planets. Therefore, the most fundamental question about these candidates is the fraction of bona fide 
planets. Estimating the rate of false positives of the overall Kepler sample is necessary to derive the 
planet occurrence rate. 

We present the results from two large observational campaigns that were conducted with the Spitzer 
Space Telescope during the the Kepler mission. These observations are dedicated to estimating the 
false positive rate (FPR) amongst the Kepler candidates. We select a sub-sample of 51 candidates, 
spanning wide ranges in stellar, orbital and planetary parameter space, and we observe their transits 
with Spitzer at 4.5 /rm. We use these observations to measures the candidate’s transit depths and in¬ 
frared magnitudes. An authentic planet produces an achromatic transit depth (neglecting the modest 
effect of limb darkening). Conversely a bandpass-dependent depth alerts us to the potential presence 
of a blending star that could be the source of the observed eclipse: a false-positive scenario. 

For most of the candidates (85%), the transit depths measured with Kepler are consistent with the 
transit depths measured with Spitzer as expected for planetary objects, while we find that the most 
discrepant measurements are due to the presence of unresolved stars that dilute the photometry. The 
Spitzer constraints on their own yield FPRs between 5-40%, depending on the KOIs. By considering 
the population of the Kepler field stars, and by combining follow-up observations (imaging) when 
available, we find that the overall FPR of our sample is low. The measured upper limit on the FPR of 
our sample is 8.8% at a confidence level of 3cr. This observational result, which uses the achromatic 
property of planetary transit signals that is not investigated by the Kepler observations, provides an 
independent indication that Kepler's false positive rate is low. 

Subject headings: planetary systems — transits — techniques: photometry 


1. INTRODUCTION 

Launched in March of 2009, NASA’s Kepler mission 
is a space-based photometric telescope designed to ad¬ 
dress important questions on the frequency and char¬ 
acteristics of planetary systems around Sun-like stars, 
and to search for transiting Earth-analogs (iBorucki et al.l 
l20I0aj) . Statistical answers to these questions are re¬ 
quired in order to constrain planetary formation and 
evolution scenarios. Kepler detects transiting plane- 
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tary candidates signals through continuous photometric 
monitoring of about 160 000 stars a t high photometric 
precision (e.g., IBorucki et al.l [Ml 0 alibi : iKoch et all 120101 : 
iGilliland et al.ll2010all . This unprecedented sample of po¬ 
tential exoplanets has become an immense resource for 
statistical studies of the properties and d istributions of 
planets around inain-seauence stars ('e-pr.lYoudinl 120111 : 
iTremaine fc Dori3l2012l: iWu fc Lithwickll2013ll . This en¬ 
semble of candidates is also necessar y for determining the 
occurrence rate of exoplanets ('e.g.. iHoward et aLll2012l : 
iFressin et al.lIMT^ . and more specifically of Earth-size 
planets in the habitable zone of their parent stars (e.g., 
Cataiizarite fc Shad 120111: iTraubI 120121 : iPetigura et all 
20131: iDressing fc Charbonneaull201^ 

The mission has led to the detection of 2740 plan- 
etary candidates during the first two year s of opera¬ 
tion (iBatalha et al.l [2^5 iBurke et al.l 1201411 . However, 
only a small subset of these candidates have been con¬ 
firmed as true planets. This is because asserting the 
planetary nature of a transit signal requires significant 
observational follow-up and computational efforts that 
are unachievable in a practical sense for every detected 
candidate. We do not expect all the signals to be due to 
planets: many astrophysical phenomena can repr oduce a 
simil ar lightcurve to that of a transiting planet (iBrownI 
l2003f) . Indeed, false-positive contamination is one of 
the main challenges facing transit surveys such as Ke- 




































2 


Desert et al. 


pier. During the last decade, ground-based surveys ded¬ 
icated to the search of transiting planets have spent 
considerable effort in confirming the pla netary nature 
of photometrically detec t ed candidates (lAlo nso et all 
200^_lBakose£^ I2007t iCollier Cameron et al.l 1200^ 
Moutou et al. 1200^ . These suryeys haye established that 
false positiyes usually outnumber true planetary systems 
by a large factor. It has been shown that from 80 to 90% 
of the candidates are false-positiyes for the most success¬ 
ful gr ound-based exoplanet suryeys (e.g., iLatham et al.l 
120091) . For these reasons, the true false positiye rate 
of Kepler remains an actiye research area because false 
positiyes can cri tically bias estimates of planet occur¬ 
rence rates (e.g.. iMorton fc Johnsonll2011l: lMortonl[2012l : 
iFressin et al.11201^ . This is the subject of the current 

paper. 

The Kepler suryey poses new challenges for dy¬ 
namically confirming (using radial yelocity or transit 
timing yariation) the planetary nature of candidates. 
This is because of intrinsic characteristics of the Ke¬ 
pler target sample such as the large number of can¬ 
didates, the candidates’ small size (presumably of low 
mass), and the faintness of the host stars. Conse¬ 
quently, we must deyelop new methods to determine 
the origin of Kepler detectable signals. One method 
consists of in-depth statistical yalidation of candidates 
by ruling o ut false positiye scenarios one-by-one (e.g ., 
BLENDER, iTorres et al.l[20nl IMTI: IFressin et all[2nTll) ; 
it fully exploit s the information from the shap e of a tran¬ 
sit lightcurye (jSeager fc Mallen-Ornelasl[200,‘ij) . The goal 
of this method is to demonstrate statistically that a tran¬ 
sit signal is more likely to be of planetary origin than to 
be a false positiye. In the case of Kepler, this is gen¬ 
erally made possible using follow-up obseryations, such 
as spectroscopy, imaging, and multi-wayelength transit 
photometry (including with Spitzer). This was demon¬ 
strated iri_Bie_ca^oRBie first yalidation of a Super- 
Earth ([Torres et al.ll20ll]) . Howeyer, each candidate yal- 
idated by this method requires intense obseryational and 
computational follow-up work. In particular, the follow¬ 
up strat egies adopted b y the Kepler team are summa¬ 
rized in iBatalha et al.l (|20I0ari , and often require sub¬ 
stantial efforts and resources. Therefore, it is impractical 
at present to apply the BLENDER method to each in- 
diyidual Kepler transit signal. Yet at the same time, 
we require the fractional yalues of bona fide planets, 
or of astrophysical false positiyes, to accurately deter- 
mine the occurrenc e of planetary systems from Kepler 
([Fressin et al.l[2M^ . 

There are currently seyeral approaches to estimate 
the False Posi t ive R ate (FPR) of the Kepler sample. 
iCoughlin et al.l ()2ni4[ ) studied the effect of contamina¬ 
tion on the FPR due to the design of Kepler itself, 
such as direct PRF (pixel response function), antipo¬ 
dal reflections, CCD cross-talks, or columns anoma¬ 
lies. The contamination sources are eclipsing bina¬ 
ries, yariable stars, and other transiting planets and re¬ 
sults in a sign ificant number of the known KOIs to be 
false positiyes. ICoughlin et al.l ()20I4[ ) performed period- 
ephemeris matching among all transiting planet, eclips¬ 
ing binary, and yariable star sources. They examined 
the full KOI list and found that 12% of KOIs are 
false-positiyes due to contamination. Other approaches 
use generic arguments about the Kepler signals to in¬ 


fer t h e oyeral l Kep l er FPR (e.g. , IMo rton fc JohnsonI 
120111: iMortonI 1201^ IFressin et al.l I2013D . There are 
also parallel attempts to estimate the FPR of tar¬ 
geted specific samples of Kepler Objects of Interest 
(KOIs). For example, studies ha ye focused on close-in 
gas giant planet candidates ('e.g.. iSanterne et'm 120121 : 
iColdn et ai.ll2012D. o r on th e multiple planet system can- 
didates (iRowe et^ I2014D . The lat ter sample contains 
less than a percent (pf false positiyes (ILatham et al.l[2MTI : 
iLissauer et al.ll20T^ 120141 ). Other methods use a proxy 
of the host stars’ mean density t o estimate the Kepler 
FPR ('e.g.. iSliski fc KiDDingll201^ . 

In this paper, we conduct two campaigns to mea¬ 
sure transit depths of KOIs with Spitzer, and combine 
these obseryations to followup studies, in order to as¬ 
sess the oyerall FPR of these samples. We adopt an 
approach that expands significantly the number of KOIs 
that are examined using multi-wayelength photometry. 
Our project focuses primarily on smaller size candidates, 
such as mini-Neptune and Super-Earth size objects, com¬ 
pared to preyious targeted sample studies. We select 
a sample of 51 candidates, measure their transit depth 
at 4.5 microns with IRAC, and combine these obserya¬ 
tions with complementary follow-up studies and infor¬ 
mation from Kepler in order to deriye the false positiye 
probability (FPP) for each object. Our method is based 
on the fact that the relatiye depth of a planetary tran¬ 
sit is achromatic (neglecting the modest effect of limb- 
darkening), but not for a blend. In contrast, a blend 
containing a false positiye, for instance an eclipsing bi¬ 
nary, can yield a depth that can yary significantly with 
the instrument bandpass and stellar temperatures. The 
amplitude of this effect increases correspondingly as the 
difference in wayelength between the two bandpasses in¬ 
creases. Since Kepler obseryes through a broad bandpass 
at yisible wayelengths, large color-dependent effects for 
false positiyes caused by the presence of blended cool 
stars can be reyealed at infrared wayelengths. We first 
applied this method by comb ining Spitzer data fr om this 
program and Kepler data in IFressin et al.l (j20II[) . 

The two Science Exploration Spitzer programs, which 
form the core of the data presented in the current 
paper, haye been an actiye part of the attempts 
to validate KOIs. About 20% of the total amount 
of Spitzer telescope time allocated for this project 
has already been used in publications dedicated to 


the confirmation or the v alidation of 22 Kep ler plan¬ 
ets. These are Kep l er-1 O c (IFressin et al 1120111). Keoler- 
11b (ILissauer et al.ll201lD. Keoler-ldb ( Buchhave et al.l 
l2011f). Kepler-18b . c (iCochran et al.l 12011), Kepler- 
19b feallard et all (120111). Kepler-20 b c Gautier et al] 
()2012D . Kepler-22b (iBorucki et ^ 1201 2D. Kepler- 
25b.c (iSteffen et al.l I2012D. Kepler-26c (iSteffen et al.l 
I2012D Kepler-32b (iFabrvckv et al.l 120)^), Kepler- 
37b (iBarclav et ^ I2013D. KeDler-49b,c (ISteffen et al.] 
I2013D. Keoler-Olb ( Ballard et al.l l20i3D. KeDler-62 e 
( Borucki et al.ll2013D. Keoler-GSb dGilliland et al.l[20T^ . 


Kepler-410A b (iVan Evlen et al.l I2014D . Kepler-93b 
([Ballard et al.ll2014) (see Table [3] for the correspondance 
between Kepler names and KOI numbers). Furthermore, 
some of the KOIs of the current study are already con¬ 
firmed or validated as planets withou t using the Spitzer 
data. In particular. IRowe et al.l (I2014D validate 851 plan¬ 
ets in multiple-planet system candidates (including 11 
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KOIs used in our study) by applying statistical ar gu- 
ments (jLatham et alJ[201lllLissauer et ani2012L[2014fl on 
the Q1-Q8 Kepler data. While the radial velocity tech¬ 
nique allowe d the confirmation of the planet ary nature of 
Kepl er-89d (jWeiss et al.ll2013D. Kepler-94 b (jMarcv et al.l 
l2014fl . Kepler-102d.e ([Marcv et ^I2ni4[l . Nevertheless, 
the present work disregard previous validation or confir¬ 
mation of individual object in order to treat the whole 
KOI list followed with Spitzer as a statistical ensemble; 
this is necessary to estimate the FPR of this sample. 
Finally, 150 hours (11%) of time from these two Ex¬ 
ploration Science Programs were used to study the at¬ 
mospheres of Kepler-detected hot Jupiters detected by 
monitoring their secondary eclipses (iDesert et al.ll2011bt 
iFortnev et al.ll201lt IDesert et al.ll2011^ 

This paper is organized as follows: we first describe 
the different types of astrophysical false positives that 
we are concerned with (Section [2]) . We then present the 
sample of candidates that were selected to conduct this 
study (Section [3]) . The Spitzer observations and results 
are presented in Sections |4] and [5l The combination of 
various observational constraints (Section [6|) allows us 
to estimate the FPR of our Kepler sample (Section [7]). 
We finally discuss the implications of our findings in Sec¬ 
tion [SJ in particular, in the context of other studies. 

2. ASTROPHYSICAL FALSE POSITIVES IN THE KEPLER 
SIGNALS 

There are a variety of astrophysical phenomena that 
can mimic the signal of a transiting planet passing in 
front of a main sequence star targeted by Kepler. These 
events are produced by additional stars falling within 
the same aperture as the target star (presumed to be 
the brighter star), and significantly diluting the total 
light observed by Kepler. More specifically, the kinds 
of false positives that we are concerned with in our 
study include background or foreground eclipsing bina¬ 
ries (EBs), blended within the Kepler aperture of the 
target star, as well as those that are physically associ¬ 
ated to the target star. We refer to these latter config¬ 
urations as hierarchical triples (HTs). HTs often cannot 
be resolved in high-angular resolution imaging. In this 
paper, we check for the presence of a stellar companion 
by looking at how transit depths vary between the Ke¬ 
pler and Spitzer bandpasses. The dilution by a stellar 
companion, blended in the Kepler aperture of a planet 
host star of interest, can be responsible for variations in 
the wavelength-dependent transit depths measured for a 
planet. We do not consider the case where the contam¬ 
inating star is itself transited by a planet as a potential 
false positive scenario. 

Before searching for false-positives in Kepler, it is 
important to recall the major vetting steps that each 
Kepler target goes through. First, a comprehensive 
study was applied when assembling the Kepler In¬ 
put C at alog (KIC: iLatham et al.l l2005t iBatalha et al.l 
l2010allbl : iBrown et al.l 1201 ili ~ leading to the identifica¬ 
tion of some EBs and stellar giants, hence avoiding 
their continuous monitoring with Kepler. About 160 000 
stars were carefully selected from the KIC catalog and 
were continuously monito re d photometrica ll y with Ke¬ 
pler (iJenkins et al.l 1201011 . IBatalha et al.l dboiObH ex¬ 
plains the detection of transit events and the vetting pro¬ 
cesses that are then applied to reject the most common 


false positive scenarios. Transit-like signals are identifi¬ 
able from TCEs (Threshold-Crossing Events) using the 
Kepler photometry alone. The Kepler team adopted a 
detection threshold of 7.1 a for the transit so that no 
more than one spurious signal can occur from purely ran¬ 
dom fluctuations amongst the 160000 stars. In practice, 
the process of vetting from TCEs to KOIs involves sev- 
eral qualitative steps tha t could affect the Kepler EPR. 
iChristiansen et aP (1201311 have validated th e integrity of 
this threshold while iCoughlin et al.l (|2014l l report that 
TCERT (Threshold Crossing Event Review Team) is 
92.9% effective in detecting false-positives for KOIs from 
Q1-Q8. 

The Kepler pipeline identifies grazing EBs by search¬ 
ing for even/odd transit depth differences or by look¬ 
ing for the presence of a clear signature of secondary 
eclipses. Giant star-eclipsed-by-a-dwarf star scenarios 
are detected by recognizing that the primary star is a 
giant, thereby implying that the size of the transiting 
body must itself be stellar (e.g., iGilliland et al.l I^OlOl : 

I Huber et al.ll2013ll . The detection of the shift in the pho¬ 
tocenter at a significance level greater than 3tT and the 
comparison of the difference of in- and out-of transit im¬ 
ages show the true source location. Interestingly, this 
technique permits the identification of potential contam- 
inati on by unresolved close-by EBs in an effi cient man¬ 
ner (| Jenkins et al.l[2M^ iBrvson et al.l 1201311 . However, 
even for high transit SNR candidates, some blended bi¬ 
nary scenarios remain undetectable through the vetting 
processes. Therefore, estimating the FPR from Kepler 
requires knowledge of the probability of encountering 
such blend scenarios. 

Throu ghout this paper, w e follow the notation intro¬ 
duced by iTorres et al.l (j20lHl : the objects that comprise 
a blended binary system are referred to as the ‘secondary’ 
and ‘tertiary’, and the candidate star host is referred 
to as the ‘primary’. The distance along the line-of- 
sight between the binary system and the main star is 
parametrized in terms of the difference in distance mod¬ 
ulus, /r. The notation applies to every astrophysical false 
positive scenario. 

3. SELECTION OF THE KEPLER OBJECT OF INTERESTS 

The first Spitzer follow-up program comprised 36 
of t he first 400 KOIs i dentified by the Kepler sur¬ 
vey (|Borucki et al.l 1201111 . A second set of 23 can- 
di dates was selected from the 2335 KOIs compiled 
bv IBatalha et al.l (|2013[1 for the second Spitzer program. 
We present the observations and the results for the 36 
KOIs from the first program and 15 KOIs from the 
second program; the current project uses in total 51 
KOIs. The remaining 8 targets from the second program 
were not observed yet at the time of the present anal¬ 
ysis (KOIs-248.03, 1686.01, 2290.01, 2124.01, 2311.01, 
2418.01, 2474.01, 2650.01). KOIs-248.03 and 2650.01 
have recently been validated as Kepler-49d, and Kepler - 
395c, respectively (iSteffen et al.ll20i^lRowe et al.ll2014ll . 
Importantly, we select these two ensembles with very dif¬ 
ferent criteria. For the first set, our goal was to derive 
the FPR of a sub-sample by following-up candidates that 
represent the diversity of KOIs initially found. To reach 
this goal, the first sample is chosen to cover represen¬ 
tative ranges of orbital periods, transit depths, stellar 
types and magnitudes that the first 400 KOIs could allow. 
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The size of this sample includes about 10% of the known 
KOIs at the time. The ranges of the second set are more 
tightly constrained: we select candidates for which the 
expected planetary radius Rp would be less than 1.6 
and the Teq < 350 K (using a stellar temperature T^e 
estimated from the KIC). Overall, our sample of stellar 
hosts span a range of KIC estimated temperature from 
Teff = 3700 K up to Teff = 9000 K. For both ensembles, 
the requirement was imposed to accept targets with pre¬ 
dicted transit depths detection of at least 3 a (scaled from 
the Kepler value), achievable with three or less transits 
observed with Spitzer. The selected candidate radii as 
function of their periods are presented in Figure [T] 

4. FOLLOW-UP OBSERVATIONS OF SELECTED Kepler 

CANDIDATES WITH Spitzer 

4.1. Spitzer observations 

We use Warm -Syitzer/TRAC (jWerner et al.l 12004 
iFazio et "ni l2004fl at 4.5 /im to observe transits of the 
51 selected KOIs between May 2010 and July 2012. We 
obtained these observations as part of two large Science 
Exploration Programs (program ID 60028 and 80117). In 
total, 1400 hours of Spitzer time is used for the follow-up 
of Kepler targets. Of this time, 800 hours of observa¬ 
tions are used to complete the first program (60028) and 
600 hours are dedicated to the second program (80117). 
150 hours of this time were used to study the atmo¬ 
spheres of hot- Jupiters detected by Kepler during sec¬ 
ondary eclipses (iDesert et aI1l2011bl : iFortnev et al.ll201ll : 
iDesert et al.ll2011cll . The remaining 1250 hours are ded¬ 

icated to validating KOIs and estimating the FPR in the 
Kepler data, and are the focus of this paper. A total of 
157 Spitzer AORs (Astronomical Observation Requests) 
have been submitted for these two programs. This paper 
focuses on the study of the 51 KOIs presented in Sec¬ 
tion [3] that have been observed with Spitzer during 95 
visits (AORs). 

For most of the targeted stars, the data were obtained 
in a continuous staring full array mode (256 x 256 pix¬ 
els) with exposure times of 12 or 30 s, depending on the 
brightness of the star of interest. We used the subarray 
mode of IRAC for the brightest host stars. In this mode, 
only a 32 X 32-pixel part of the detector is used; this 
covers a 38 x 38 arcsec^ field of view (pixel size of 1.2") 
and allows for higher cadences (0.2 s exposures). We 
choose to put our target at the default pointing position 
in the center of the field-of-view in order to avoid known 
hot pixels and bad columns. This area of the detector is 
well characterized since it has been extensively used for 
extrasolar planet studies. An offset is applied to a few 
KOIs in order to avoid the presence of a close-by bright 
target on the same line or column. The ephemerides 
of the KOIs were taken from the KFOP database, which 
are now available on the CFOP websit43 and we ensured 
that each visit lasted approximately 2.5 times the transit 
duration. We observe 29 KOIs, amongst the 51 presented 
here, during multiple transit events in order to improve 
the SNR on the combined lightcurves. Tables [T] and [5] 
list these observations for each program respectively. 

4.2. Spitzer photometry 

https://cfop.ipac.caltech.edu/home/ 


We use the BCD files (Basic Calibrated Data) pro¬ 
duced by the Spitzer/IRAC pipeline. These files in¬ 
clude corrections for dark current, flat fi elding, pixel 
non-li nearity, and conversion to flux units. iDAsert et al.l 
(l2009fl describes the method used to produce photomet¬ 
ric time series in each channel from the BCD files. The 
method consists of finding the centroid position of the 
stellar point-spread-function (PSF) and performing aper¬ 
ture photometry. We first convert the pixel intensities to 
electrons based on the detector gain and exposure time 
provided in the FITS headers. This facilitates the eval¬ 
uation of the photometric errors. We extract the BJD 
date for each image from the FITS headers and com¬ 
pute it to mid-exposure. We then correct for transient 
pixels in each individual image using a 20-point sliding 
median filter of the pixel intensity versus time. To do 
so, we compare each pixel’s intensity to the median of 
the 10 preceding and 10 following exposures at the same 
pixel position and we replace outliers greater than 4 a 
with their median value. The fraction of pixels that we 
replace varies between 0.15 to 0.5%. 

The centroid position of the stellar PSF is then deter¬ 
mined using DAOPHOT-type Photometry Procedures, 
GCNTRD, from the IDL Astronomy LibrarjQ- We use 
the APER routine to perform aperture photometry with 
a circular aperture of variable radius. For each visit, we 
search for the best aperture size ranging between 1 and 
8 pixels radii in steps of 0.5 pixel. We propagate the 
uncertainties as a function of the aperture radius and 
we adopt the size that provides the smallest errors. We 
notice that the SNR does not vary significantly with the 
aperture radii for all the dataset. The final aperture sizes 
are set between 2.5 and 3.5 pixels depending on the KOIs. 

We determine the background level for each frame from 
two methods. The first method uses a fit of a Gaussian to 
the central region of a histogram of counts from the full 
array, where the background values are defined by the 
peak position of this Gaussian. The second method uses 
the measure of the median value of the pixels inside an 
annulus centered around the star, with inner and outer 
radii of 12 and 20 pixels respectively, to estimate the 
background overall level. Both estimates produce similar 
results. The contribution of the background to the total 
flux from the stars is low for all observations, from 0.1% 
to 1.0% depending on the images, and fairly constant for 
each AOR. We find that the residuals from the final light 
curve modeling are minimized by adopting the center 
of the Gaussian fits. After producing the photometric 
time-series, we use a sliding median filter to select and 
trim outliers greater than 5cr, which correspond to less 
than two % of the data. We also discard the first half- 
hour of all observations, which is affected by a significant 
telescope jitter before stabilization. 

Six AORs, corresponding to a total of 30 hours, were 
gathered at levels above 30,000 DN (in the raw data); 
this is a level where the detector response tends to be 
non-linear by several percent. In order to avoid misin¬ 
terpreting these data, we do not consider AORs that are 
above the range of linearization correction. Therefore, 
these six AORs are not used in this work and they are 
not presented in Table [I] and [5] 
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4.3. Determination of the transit depths from Spitzer 
lightcurves 

As described in iDesert et al.l (j2011a[l , we use a tran¬ 
sit light curve model multiplied by instrumental decor¬ 
relation functions to measure the transit parameters and 
their uncertainties from the Spitzer data. We compute 
the transit light c urves with the I DL t ransit routine 
DCCULTSMALL from iMandel fc Agoll (120021 1 . This model 
depends on the following parameters: the planet-to-star 
radius ratio the orbital semi-major axis to stellar- 

radius ratio (system scale) a/i?*, the impact parameter 
b, the time of mid transit Tc, and limb darkening coeffi¬ 
cients. 

The measured parameter of interest here is the transit 
depth. We hx Tc, a/Ri, and b to their values measured 
from the Kepler photometry. The SNR of our observa¬ 
tions is low compared to typical Spitzer observations of 
brighter transiting planets. The limb darkening effect is 
negligible at this level of precision; the coefficients are 
set to zero. We assume an eccentricity of zero for the 
KOIs orbit since this parameter does not affect the tran¬ 
sit depth measurements at the level of precision we are 
working with. Only is set as a free parameter to 

represent the astrophysical signal. 

The Spitzer/IRAQ, photometry is known to be system- 
atic ally affected by the so-cal led pixel-phase effect (see 
e.g.. iCharbonneau et al.l[2?i0^ . This effect is seen as os¬ 
cillations in the measured fluxes with a period of ap¬ 
proximately 70 minutes (period of the telescope pointing 
jitter) for data secured prior October 2010, and 40 min¬ 
utes for data secured after. By October 2010 the Spitzer 
engineering team was able to correlate the pointing wob¬ 
ble with the cycling of a heater used to keep a battery 
within its operating temperature range. Following this 
discovery, the Spitzer team signihcantly reduced the am¬ 
plitude and the period of the pointing wobble. The am¬ 
plitude of this effect varies between 1 and 2%, peak to 
peak, depending on the position of the star in the ar¬ 
ray. We decorrelated our signal in each channel using a 
linear function of time for the baseline (two parameters) 
and three types of functions to correct the data for the 
intrapixel variations: a linear function of the PSF po¬ 
sition (two parameters), a quadratic function (four pa¬ 
rameters) a nd a quadrat i c wit h a cross term (five pa¬ 
rameters). IDesert et al.l (|20n9[ l describes in detail the 
last function. 

We perform a simultaneous Levenberg-Marquardt 
least-squares fit (|Markwardtll2009ll to the data to deter¬ 
mine the transit depths and the instrumental parameters. 
For each visit, we adopt the decorrelation function that 
significantly improves the minimization. We rescale 
the errors on each photometric point to be set to the root- 
mean-square (rms) of the residuals from the initial best- 
fit of the data. Hence, the reduced becomes one. All 
the data-point measurement errors are therefore assumed 
to be identical for each lightcurve. As an example. Fig¬ 
ure [2] shows the raw data and the corrected Spitzer tran¬ 
sit lightcurve of KOI-701.03 (Kepler-62e). Figures |3l HI 
and [5] present the normalized, corrected, binned and 
combined lightcurves with their associated best fit mod¬ 
els for all the observed KOIs that are presented in the 
current study. 

We estimate parameter uncertainties using two differ¬ 


ent methods: Markov Chain Monte-Carlo (MCMC) and 
residual permutation methods. Our MCMC implemen¬ 
tation use s the Metropolis-Hasting algorith m with Gibbs 
sampling (jTegmark et al.l[?nn4 Ih^ordll2005fl . We assume 
uniform prior distributions for all jump parameters. We 
adjust the width of the distribution from which we ran¬ 
domly draw the jump sizes in each parameter until 20- 
25% of jumps are executed in each of the parameters. 
We create five chains, each with 10^ points, where each 
chain starts with a different set of starting parameters 
(each parameter is assigned a starting position that is 
+3 ct or —3tT from the best-fit values). We discard the 
first 10% of jumps of each chain to remove the chain’s 
transient dependence on the starting parameters. 

In order to obtain an estimate of the correlated and 
systematic errors in our measurements, we use the resid¬ 
ual permutation boots trap, also ca ll ed “Pr ayer Bead” 
method, described in iDasert et al.l (l2011ari . In this 
method, the residuals of the initial fit are shifted system¬ 
atically and sequentially by one frame, and then added 
to the transit light curve model before fitting again. 

For both methods, the posterior distributions are used 
to estimate the errors: we allow asymmetric error bars 
spanning 34% of the nearest points above and below the 
values of the parameters associated with the minimum x^ 
to derive the 1 a uncertainties for each parameter. We 
find that the two approaches provide consistent results. 
Tables [T] and [2] present the transit depths and associated 
errors derived from the MCMC technique. 

Finally, we check that KOIs for which we have multi¬ 
epoch measurements have transit depths that agree 
within the 3 a level. We combine the measured transit 
depths for these KOIs by computing the weighted means 
and errors. 

4.4. Determination of the Spitzer magnitudes 

We use standard aperture photometry of each individ¬ 
ual BCD image in order to compute the flux for all the 
KOIs in our sample. We measure the averaged flux over 
the background annulus. The main difference between 
the proced ure used in the section and the one described 
in Section 14.31 is that we use a fixed aperture size with 
a radius of 3 pixels surrounded by an annulus of I2-to- 
20 pixels to estimate the flux and the background, re¬ 
spectively. Furthermore, only the out-of-transit data are 
considered for determining the source flux densities. We 
first convert the BCD images into mJy per pixel units 
from their original MJy per steradian units. Then we 
estimate the centroid position of the main target star’s 
PSF for each image, and compute aperture photometry 
centered on the source. We apply an aperture correc¬ 
tion of 1.113 at 4.5/rm (this value is taken from IRAC 
data handboolO). We correct the full lightcurve for the 
intrapixel sensitivity using the method described in Sec¬ 
tion 14.31 Color and array-location-dependent photomet¬ 
ric corrections are also applied to the photometry; the 
latter accounts for the variation in pixel solid angle (due 
to distortion) and the variation of the spectral response 
(due t o the tilted filters and wide field-of-view) over the 
array (iHora et al.ll2004ll . Most of our data taken in full- 
array mode is such that the PSF is centered on the cen¬ 
tral pixel of the array (128;128), so no array correction 

15 
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is applied for this dataset. The brightest stars of our 
sample are secured in subarray mode for which we ap¬ 
ply an array correction of 0.68%. We convert the surface 
brightness (in mJy) into Spitzer magnitudes at 4.5/xm 
using a zero-m agnitude flux densit y (zmag) of 179.7 Jy 
as computed bv lR.each et al.l ()2nn,^ . We finally compute 
the uncertainties on the flux densities using the photon 
noise and the standard deviation of the measurements 
from all the individual frames. We test the accuracy 
of our procedure using Spitzer IRAC photometric cal¬ 
ibrator datasets taken from the public Spitzer archive: 
BD-l-60 1753 and HD 180609. We c heck t hat our magni¬ 
tudes match those of iReach et al.l (|2005f) at better than 
the 1 a level. Tables [T] and [ 2 ] present the flux densi¬ 
ties and the corresponding magnitude at 4.5 /xm for each 
KOI. While the uncertainties in these tables are the for¬ 
mal values, we have conservatively adjusted some of the 
errors to be no lo wer than 2%. This l ower limit is based 
on the findings of IReach et al.l (1200511 . 

5. USING Spitzer OBSERVATIONS TO RULE-OUT FALSE 
POSITIVE SCENARIOS 

5.1. Analytical framework 

This section describes in detail how we use Spitzer ob¬ 
servations to rule-out false positive scenarios. The ap¬ 
plied methodology makes use of the transit depths mea¬ 
sured with Spitzer and with Kepler, as well as the mea¬ 
sured colors Kepler-Spitzer. 

The true transit depth St obtained from an eclipsing 
system comprising a main object ( 2 ) and its companion 
(3) corresponds to: 




SF2 

F 2 + F 3 ’ 


( 1 ) 


where F 2 and F^ are the emitted fluxes in the same band¬ 
pass. The parameter S represents the surface ratio be¬ 
tween the two objects and is expressed as: 


S = 



( 2 ) 


The blended transit depth corresponds to the apparent 
transit depth Sb of this eclipsing system diluted with a 
primary star (1). It is computed as follows: 


5b 


5 F 2 _ 

Fi + F2 + F3 


(3) 


where Tj is the flux from the primary star, the targeted 
KOI, and d is the dilution due to the presence this star 
in the line of sight of the eclipsing binary. 

The ratio of the apparent transit depths measured in the 
Kepler (K) and the Spitzer (S) bandpasses corresponds 
to: 


5b,s _ 5t,s ds . , 

5b,K 5t,K dK ’ 

The ratio of the true transit depths is calculated from 
Eq.m and is expressed as: 


In order to simplify the problem, we assume that the 
contribution of the tertiary flux to the ratio of the dilu¬ 
tion in the two bandpasses is negligible. Omitting the 
contribution of the tertiary has a similar effect as reduc¬ 
ing the distance modulus. Therefore, this approximation 
does not significantly impact the final results. Under this 
assumption, the ratio of the dilution in the Spitzer band¬ 
pass to that in the Kepler bandpass can then be written 
as: 


dj^ 10-0-4(Mi,s-M2,s-|-Ai) -(- 1 ’ 

where Mi and M 2 are the absolute magnitudes of the 
primary and the secondary stars. The distance between 
the binary and the primary star is parametrized for con¬ 
venience in terms of the difference in distance modulus, 
/X (/X = 0 for HT scenarios). 

Finally, the three component “Kepler—SpitzeF’ color 
C'kSi 23 i® expressed as: 

C'kSi23 =-2.5 X 

10~O.4(Mi,k-I-m) _|_ J^Q—0.4(M2,k) _j_ 0.4(M3,k) 

_|_ ]^Q —0.4(M2,s) _|_ ;[Q—0.4(M3,s) 

(7) 


5.2. Using the analytical framework to explore blend 

effects 

As a first step, we apply the above analytical frame¬ 
work to reveal the effect of blends in various scenarios 
that we expect to encounter. In particular, we explore 
blend types involving EB and HT systems that can mimic 
the KOI properties in both Kepler and Spitzer band- 
passes. 

We simulate blends us ing model isochrone s from the 
Padova isochrone series (iGirardi et al.l [200^ . This al¬ 
lows us to set the properties of the three stars involved, 
specifically their masses, from which their brightnesses 
can be predicted in any bandpass. The relevant band- 
passes for this work are those that correspond to Kepler 
and Spitzer. These bands are available on the Padova 
models’s CMD websitfF^. For simplicity, we adopt here 
a representative isochrone of 3 Gyr and solar metallicity. 
The validity of this approximation, which has only a mi¬ 
nor impact on our final results, is discussed in Sect 18. 3 1 

We now present some of the characteristic features 
from our simulated blends, and we discuss how they 
change with the adopted stellar parameters. Figure [ 6 ] 
displays the ratio of the true eclipse depths in the Spitzer 
and Kepler bands for an undiluted eclipsing binary 
(eq.m as a function of the mass of the tertiary star (M 3 ). 
We show this ratio for three different masses of the sec¬ 
ondary star (M 2 ). The net dilution effect caused by the 
primary star as a function of secondary mass for three 
values of Mi is presented in Figure!?] Figure |S| shows the 
same three cases as in Figure IH though we include this 
time the effect of the dilution in the signals produced by 
a primary star of one solar mass (Mi = IMq). These 
apparent Spitzer/ Kepler eclipse depth ratios (eq. [ 6 ]) are 
shown for two different blend scenarios: a HT configura¬ 
tion (in which the difference in distance modulus between 


5t,s _ _ F 2 ,s/F 2 ,K _ 

St,K (7^2,3 + Fa^s)/(F 2 ,k + Fs^k) 


(5) 


4® http://stev.oapd.inaf.it/cgi-bin/cmd 
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the primary and the EB is zero), and a configuration with 
the EB in the background (a distance modulus of five). 
Finally, figure [9] illustrates how the stellar properties af¬ 
fect the color difference of a blended system for colors 
computed between the Kepler and Spitzer bandpasses. 
This plot shows the difference between the color of the 
combined light of the primary and secondary {C 12 ) and 
the color of the primary alone (Ci), as a function of the 
secondary’s mass. For this illustration, we assume that 
the tertiary’s contribution is negligible, as will usually be 
the case. 

5.3. Applying the framework to the Spitzer and Kepler 
observations 

The methodology described above allows us to reject 
many false positive scenarios for the KOIs observed in 
this study. The observational constraints consist of the 
measured transit depths and the apparent magnitudes 
in the Spitzer and Kepler bandpasses (see Tables 1 and 
2). The tests for potential blends use eqs. 5 and 6 . The 
free parameters considered are the secondary and tertiary 
masses, and the relative distance between the eclipsing 
pair and the primary star (i.e., the difference in distance 
modulus). We explore these quantities over wide ranges 
of stellar masses and distance moduli, in a grid pattern 
to fully map the space of parameters for allowed blends. 
Primar y masses for th e KOI s (Mi) are taken from the 
work of iBatalha et al.l (l2013f) and held fixed. We allow 
M 2 and M 3 to vary between 0.1 Mq and 1.4 Mq, and the 
distance modulus difference p is divided in linear steps 
between —5 mag and 15 mag. For each star, the intrinsic 
brightness in the Spitzer and Kepler bands is read off di¬ 
rectly from the adopted isochrone. At each trial distance 
modulus we compute the ratio between the transit depth 
from Spitzer (- 1 - Scr) and that from Kepler. This ratio 
sets an upper limit to the mass of the secondary that 
could mimic the transit depths measured in both ban- 
passes. We do not set a lower limit between the transit 
depths from Spitzer (— 3cr) and that from Kepler, be¬ 
cause these shallower depths would involve scenarios with 
massive stars (these are typically more massive than the 
Sun; See Figure [H). In practice, the number of stars that 
we eliminate from the shallower limits represent only a 
small fraction of the number of low mass stars removed 
from the deeper limits. In a typical Kepler aperture, the 
Besangon population synthesis model shows that only 
10% of the stars have mass greater than 80% of the Sun. 
Consequently, we only consider the deeper transit depth 
from Spitzer (- 1 - Scr) in our calculation of the FPPs. 

Similarly, we use the Kepler minus Spitzer color to 
further constrain the blend properties. The color is de¬ 
rived from our measured Spitzer magnitude (1115 ± Scr 
; Table [J and and from the Kepler magnitude as re¬ 
ported bv iBatalha et al.l (I2013D . Assuming that the mea¬ 
sured fluxes result from the contribution of three stars, 
the color €123 (eq. [7|) is then compared with the color of 
the primary star alone, Ci, computed from the adopted 
isochrone. Only a subset of the secondary and tertiary 
stars can reproduce this color difference, which provides 
another constraint on the secondaries’ masses. 

6. ADDITIONAL OBSERVATIONAL CONSTRAINTS 

This section presents additional observational con¬ 
straints used for several candidates of our sample. These 


complementary observations allow us to exclude more 
false positive scenarios, which remain after applying the 
constraints from the Spitzer observations. 

6 . 1 . Stellar Reconnaissance 

In general, follow-up observations of Kepler planet can¬ 
didates involve reconnaissance spectroscopy. This is nec¬ 
essary to characterize the primary star a nd to look for 
eviden ce of astrophysical false positives (jBatalha et al.l 
I2010al) . Such false positives include single- and double- 
lined binaries, some HTs and EBs, which would show ve¬ 
locity variations at large amplitudes. We also use their 
spectra to estimate the effective temperature, surface 
gravity, metallicity, and rotational and radial velocities of 
the host star. For the current study, we assume that the 
primary star is the brightest star and that it is known and 
characterized. In theory, if we are to consider all possi¬ 
ble scenarios, a secondary star could be brighter than the 
primary target (see Figure [5]). In practice, our assump¬ 
tion that the primary star is characterized means that 
the secondary stars considered in this study can only be 
fainter than the primary stars. We treat the stars in our 
sample in a uniform ma nner as we rely o n the stellar radii 
and masses pr ovided in IBatalha et all (1201311 . which are 
from the KIC (jBrown et al.ll2011ll . This is a critical step 
to fix the primaries’ stellar characteristics in our FPP 
calculations. 

6 . 2 . Imaging 

Kepler's photometric aperture is typically a few pixels 
across with a scale of 3'.'98 per pixel. Therefore, high- 
resolution imaging is often performed in order to identify 
neighboring stars that may be blended EBs contaminat¬ 
ing the primary target photometry. Only 23 amongst the 
51 stars of our study have high spatial-resolution adap¬ 
tive optics images. These images were t aken in the near- 
infrar ed (J and K bands) with A RIES (iMcCarthv et alj 
on the MMT and PHARO (iHavward et al.ll20nill 
on the Palomar Hale 200-inch. The obs ervations and 
their sensitivity curves are presented in lAdams et al.l 
(1201211 . The AO images allow us to detect companion 
stars as close as 0 . 1 " from the target’s primary star. 
These images also rule-out stellar companions within a 
6 " separation from the primary, with a magnitude differ¬ 
ence up to 9. There are several KOIs for which we de¬ 
tect additional stars within the Kepler apertures of the 
primaries. We note that since we completed the analy¬ 
sis of the data for this work, several projects have been 
conducted to search for close-by stellar companions that 
could be the s ources of false positives using high reso¬ 
lution imaging (fLillo-Box et alJ[2012HAdams et al1l2013t 
ILillo-Box et al.ll2014ll . Furthermore, a robotic AO surve y 

of nearly 715 KOIs was conducted bv ILaw et al.l (I2014I1 . 
to search for stellar companion in a systematic manner. 
So far these searches found 7 KOIs within our samples 
that have detected fainter close-by stellar companions 
within 5"; these are KOT12, 13, 94, 98, 111, 174, and 
555. We also note that the environment of KOT854 has 
been scrutinize d using HST/WFC3 , but no companions 
were reported (IGilliland et al.ll20T5ll . 

6.3. Centroid analysis from Kepler 

The very high astrometric precision of Kepler allows 
us to monitor the motion of the target’s photocenter. 
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This provides an effective way of identifying false posi¬ 
tives that are caused by EBs falling within the aperture. 
We directly measure the source location via difference im¬ 
ages to search for impostors based on scrutinizing pixels 
in the KOIs’ aperture. Difference image analysis is con¬ 
ducted based on the difference between average in-transit 
pixel images and average out-of-transit im ages. A fit of 
the K epler pixel response function (PRF; iBrvson et al.l 
l2010tl of both the difference and out-of-transit images di¬ 
rectly provides the transit signal’s location relative to the 
host star. The difference images are measured separately 
in each quarter. We estimate that the transit source lo¬ 
cation is the rob ust uncertainty - weigh ted average of the 
quarte rly results. iJenkins et al.l (|2010( 1 and IBrvson et al.l 
(j2013f) describe this technique^ _ 

Subsequent to our study, IBrvson et al.l (1201,111 pre¬ 
sented the centroid analyzes for the complete list of KOIs. 
However, at the time of this study, the centroid analyzes 
have been performed only for a subset of the KOIs tar¬ 
geted in this work (see Table [3]). This analysis shows no 
significant offsets during transits in any quarter, the com¬ 
puted offsets are well within the radius of confusion (at 
the 3(7 level). This shows that the observed centroid loca¬ 
tions are consistent with the transit occurring at the KOI 
locations. Stars located at distances beyond the confu¬ 
sion radius from the targeted KOI are ignored since the 
centroid analysis would be confused by such stars and 
this would not yield accurate measurements. However, 
we note that bright stars can have PRF wings than can 
extend to ov er 100" and could ther efore contaminate the 
main target (iCoughlin et al.l[2014[l . 

6.4. Faint limit condition 

The faint limit condition corresponds to the faintest 
blended EB system than can reproduce the transit signal. 
The blended system must comprise more than a fraction 
5b (as defined in Equation [3]) of the total flux within the 
Kepler aperture. This condition may be expressed as the 
following: 

W2iim,K - wqK = Amif = -2.5 logio(^b), (8) 

where TO 2 iim.K is the apparent Kepler magnitude of the 
blended binary system and mi^K is the magnitude of the 
Kepler targeted star. This limit is such that no binary 
system fainter than m 2 iim,K can possibly mimic a transit 
signal with a depth of 5 given the primary star of Kepler 
magnitude mi^K- 


7. FALSE-POSITIVE PROBABILITY 


In this section, we compute the false positive probabil¬ 
ity for each KOI in our selected sample. For any candi¬ 
date, the rate of false positives relative to the rate of any 
event can be written as: 


FPP = 


Ffp 

Ffp + Fp ’ 


(9) 


where Ffp is the estimated frequency of false positive 
scenarios (which depends on the local stellar density), 
and Fp is the expected frequency of transiting planets 
for a given KOI. Stars that are unable to reproduce the 
observables are removed using complementary observa¬ 
tional constraints (e.g., centroid analysis, AO imaging, 
etc...). We then compute the planet and false positive 


frequencies as described in the following sections and we 
finally derive the FPP for each object. We present the 
FPP for each KOI observed with Spitzer in Table |3l Fur¬ 
ther details about the steps undertaken are given below. 


7.1. Planet Frequency 

To estimate the likelihood of a planet, we rely on the 
list of candidates from Kepler. We assume here that all 
the KOIs are true planets and that the period distri¬ 
bution of our Spitzer sample follows, in first order, the 
period distribution of the KOIs. The latest assumption 
is ve rified in Figure ITUl for the KOI distribution derived 
from iBatalha et al.l (I20I3I1 . The KOIs are separated per 
period range using equal logarithmically spaced intervals. 
We count the number of candidates that each period bin 
contains and divide by the total number of stars followed 
by Kepler (156,453) to derive the planet frequency for a 
given bin. Figure [TT] shows the distribution of planet fre¬ 
quencies that we use in the calculation of the FPP for 
the KOIs we observed with Spitzer. 

7.2. False-Positive Frequency 

In order to derive the Ffp of a KOI, we must assess the 
likelihood of the various types of false positives, and also 
estimate the local stellar density. We estimate the stellar 
density using a stellar popula tion synthesis mod el of the 
Galaxy, the Besangon model (IRobin et al.ll2003ll . We use 
this model to derive the frequency of stars present in the 
photometric apertures around each KOI in our sample. 
We adopt the stellar densities predicted by this model 
in the R band, which is a band sufficiently close to the 
Kepler bandpass for our purposes. Instead of estimat¬ 
ing the stellar population in a cone around the line of 
sight of each KOI, we create a grid of 24 cells evenly 
spaced over Kepler's field of view (about one cell per 
Kepler CCD module). An aperture of I square degree 
centered on each cell of the grid is chosen for the initial 
estimate of stellar population. We then perform stellar 
density calculations in half-magnitude bins of apparent 
brightness, from a R-magnitude of 5 down to magnitude 
24. We also account for interstellar extinction with a 
coefficient of 0.5 mag kpc“^ in V band. The number 
of stars that fall in the aperture of each grid varies be¬ 
tween 30,000 and 1,400,000, depending on the Galactic 
latitude and longitude of the KOIs. We derive the stellar 
population using the closest cell relative to the target. 
In this way, we evaluate the expected number density of 
neighboring stars, and their mass distribution per square 
degree and per magnitude bin. We then estimate the 
number of EBs amongst these neighbor stars that could 
potentially reproduce the signal detected in Kepler. It 
corresponds to the occurrence rate of binaries multiply 
by the probability that they undergo e clipses. This rate 
has been measured in the Kepler field (iPrsa et al.ll2011l : 
iSlawson et al.ll20lTll . In particular, 1261 detached binary 
systems have been found from the first fou r months of ob¬ 
serva tions with periods less than 125 days (ISlawson et al.l 
[Ml . These are the typical false positives of interest for 
our study. Therefore, the frequency of EBs that we con¬ 
sider is /eb = 1261/156,453 = 0.80%. This frequency 
is computed from short period detached binaries and it 
depends on their period; it decreases below 0.80% for 
longer periods. Consequently, the uniform frequency of 
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EBs that we use in this work is a conservative value. 
This allows us to compute upper limits for the FPPs. 
We note that this frequency includes eclipsing pairs in 
triple or higher multiplicity stellar systems. Finally, the 
total blend frequency Fpp that is inserted in Equation |9] 
is found by combining the probabilities associated with 
all background or foreground binaries. 

7.3. Combining observational constraints to remove 
blended stars within the Kepler apertures 

Eor each KOI in our sample, we use the density of 
stars in each magnitude bin determined using the popu¬ 
lation synthesis model. We now calculate the fraction of 
these stars that would remain false positive candidates 
after taking into account the available observational con¬ 
straints. We describe below the method used in this task. 

The first step uses the faint limit condition described 
in Section 16.41 This removes all the stars that are too 
faint to reproduce the transit depths observed in the Ke¬ 
pler photometry. The remaining stars pass through the 
second step that combines high-resolution images of the 
target’s neighborhood (AO observations) and centroid 
motion analyzes when available. This provides the spa¬ 
tial extent considered for blend frequency calculations. 
These constraints significantly decrease the size of the 
apertures that could hide a background or foreground 
binary: this is sometimes 5” and often 2”. Even though 
stars located at distances beyond 5” can reproduce Ke¬ 
pler signals (|Coughlin et al.lT2ni4l l. we assume for sim¬ 
plicity that all the stars that are beyond the limits found 
by the centroid method cannot mimic transit signals. 

These stars are therefore excluded from the blend 
frequency calculations. This step always removes the 
largest number of possible blend scenarios. It is a very 
powerful tool fo r identifying backgro und EBs blended 
with the target (jBatalha et al.l[20I0all . In a final third 
step, we use the constraints from the Spitzer observa¬ 
tions, the measured transit depths and magnitudes, as 
described in Section [S] For each secondary star that 
survives the first two steps, we consider tertiary stars 
(EBs), ranging from 0.1 to 1.4 solar masses, and test the 
ratios of transit depths. We compute the diluted ratio of 
the true transit depths for each scenario following Equa¬ 
tion m We compare this ratio to its observed value for 
each mass of the chosen tertiary star. If the calculated ra¬ 
tio is not consistent with the observed one, the tertiary 
star scenario is then rejected. We also apply the con¬ 
straint from the color Kepler-Spitzer" C'kSi 23 given 
by Equation [71 We reject the tertiary star scenarios for 
which the calculated colors are not in agreement with 
the measured ones. Overall, this third step allows us to 
remove most of the red dwarfs that could potentially re¬ 
main as false positives. We finally find the total blend 
frequency by combining the probabilities associated with 
all background or foreground star. 

8. DISCUSSION 

8 .1. Comparing the Kepler and Spitzer transit depths 

Eigure [T2| (top panel) shows the measured transit 
depths in the Spitzer bandpass compared with those 
measured with Kepler. We find that 50% of the sam¬ 
ple have measured depths that agree within Icr, and that 
85% agree within 3cr (Figure [T2l bottom). The distribu¬ 


tion is therefore somewhat broader (by ^20%) than ex¬ 
pected for a Gaussian with a standard deviation of unity, 
indicating that the Spitzer and Kepler transit depths for 
our KOIs are not all statistically consistent within their 
uncertainties. This may be caused by: (i) the presence of 
false positives in our sample; (ii) dilution from unresolved 
companion stars resulting in wavelength-dependent tran¬ 
sit depths; or (iii) underestimated uncertainties in our 
measurement of the Spitzer transit depths. Below we 
present evidence that some of the KOIs do indeed suffer 
from dilution effects. Likewise, biases in determining re¬ 
liable uncertainties for our Spitzer measurements cannot 
be entirely ruled out, as they depend on our ability to 
correct the data for the main source of systematic errors, 
which is the intra-pixel sensitivity. 

There are three candidates for which the transit depths 
measured from the Kepler photometry are significantly 
deeper (> 4 cr) than the depths measured from Spitzer: 
KOI-I2.01, KOI-I3.01, KOI-94.01. These objects are 
three Jupiter-size planet candidates out of the four from 
this family that we have in our sample. Since we observed 
KOI-I3.01 and KOI-94.01 with Spitzer, the se candidates 
have b een c onfirmed as bona fi de planets bv iBarnes et al.l 
(|20lH) . and I Weiss et al.l ()20I3[) respectively. These three 
KOIs have known close-by companions which are lo¬ 
cated at closer than I”, well within the Kepler aper¬ 
tures and within the Spitzer PSFs. The host star KOI- 
13 is known to be part of a stellar binary s ystem, both 
components being rapid ly rotating A stars (jSzabo et al.l 
l20 III: [Barnes et ^l20IIfl . The companion is about 300 K 
cooler than the primary host star and is 0.3 magnitude 
fainter in the Kepler bandpass. Similarly, the CFOP 
shows that KOI-12 is also a massive fast rotator star. 
Direct images of the close environment of KOI-12 reveal 
the presence of two fainter stars within 1” of the primary 
host. A low mass-star companion at 0"6 from KOI-94 has 
been detected and this expl ains the significant diff erence 
of measured transit depths (iTakahashi et al.ir2013ll . The 
candidates with the largest discrepancies are KOI-82.01 
and KOI-98.01 at the 3.7 and 3.6cr level, respectively. 
For KOI-98.01, the ho st star is known to have a stel¬ 
lar co mpanion at 0"3 (jBuchhave et al.l 120111 : iLaw et ahi 
120141). KOI-82 has no close by companion detected by 
AO (iMarcv et al.l 120141 : we therefore attribute the dis¬ 
crepancy between Kepler and Spitzer to statistical fluc¬ 
tuations. In general, the dilution produced by the pres¬ 
ence of stars within the aperture of the KOI result in 
chromatic differences between the transit depths. In the 
present cases, the flux contamination from the compan¬ 
ion stars to the primary host targets, KOI-12, 94 and 
13, vary with wavelength. The dilution is greater in the 
infrared compared to the visible as expected for cooler 
(redder) contaminants. Because the contaminant stars 
contribute proportionally to more flux in the infrared 
compared to the visible, the dilution produced by their 
presence is larger at longer wavelengths. This results in 
shallower measured transit depths in the Spitzer band¬ 
pass compared to the Kepler bandpass. 

KOI-258.01 is the only observed candidate for which 
the lightcurve is not properly fitted by a transit planet 
model. The Spitzer and the Kepler transit lightcurves 
exhibit a clear V-shape. Such objects do not usually 
appear in the KOI list as they are flagged as false posi¬ 
tives. KOI-258.01 was ranked as a KOI early on in the 
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mission before being removed. The Kepler project has 
since changed this target to an inactive mode as the pri¬ 
mary star has two companions within 1.5” that are 2 to 
3 magnitudes fainter in the K band. However, by the 
time of this decision the order to execute the Spitzer ob¬ 
servation had already been given. For consistency we 
therefore choose to remove this object from the list of 
KOIs with computed FPPs. 

8 .2. On the False Positive Rate 

As described in Section[71 we estimate the FPP for each 
individual candidate on this list and Table [3] presents our 
results. We find that half of the targets in our sample 
have a FPP which is lower than 1% (see Fig. [T4l) . Using 
the distribution of FPP for the complete sample, we mea¬ 
sure a median value of 1.3%. We calculate a robust esti¬ 
mate of the dispersion of the FPP distribution using the 
median absolute deviation as the initial g uess, and then 
weigh ting points using Tukey’s Biweight (iHoaglin et al.l 
1198311 . The dispersion of our distribution measured by 
this method is 2.5%. This leads to an upper limit for the 
FPR of 3.8% at the Icr level, and of 8.8% at the 3cr level. 

At present, the sample of this study repr esents only 2% 
of th e total candidates published so far (jBatalha et al.l 
1201311 . Therefore, extrapolating a Kepler FPR from such 
a small sample should be done with caution. Further¬ 
more, our sample is not a uniform sample of KOIs as it 
can be divided into two categories. The first category 
comprises the sub-sample of two thirds of the targets for 
which we have complementary observational constraints 
from direct imaging (AO) or KepleFs centroid measure¬ 
ments. For this first sample, we measure a median FPR 
of 0.7% with a robust estimate of the dispersion of 0.8%. 
The second category comprises the rest of the targets for 
which we have no other constraints other than the Ke¬ 
pler depths and the Spitzer photometry. This is because 
we do not have the centroid analysis completed yet at 
the time of the work. For this second sample, we mea¬ 
sure a FPR lower than 24% at Icr level, and we find that 
these FPPs vary between 5-to-43 %. The KOIs of our 
sample with the highest FPP values (e.g. KOT137.01, 
137.02 and 952.02) are indeed confirmed planets, which 
implies that their FPP is much lower than the values we 
can compute from the current study. All of this leads to¬ 
wards a rather low value for the overall FPR of the Kepler 
sample. Finally, we note that 11 of the observed KOIs 
have FPPs < 0.3%. This implies that these detected sig¬ 
nals are at least 99.7% consistent with planetary objects; 
these KOIs can be considered as “validated” at the 3cr 
level of confidence. 

We use two Spitzer Science Exploration Programs for 
this study. We measure a median FPR of 0.8% with a 
dispersion of 0.9% from the KOIs of our first Spitzer pro¬ 
gram (60028). This value becomes 14%, with a disper¬ 
sion of 9%, for the 15 targets from our second program 
(80117). This is a higher FPR compared to the value 
derived from the first Spitzer program. This is mainly 
because we do not have constraints from AO imaging nor 
from the Kepler centroid analysis for these KOIs, un¬ 
like for most of the targets from the first program. This 
higher FPR value is also because the targeted sample is 
mainly composed of M-dwarf hosts. Low mass stars are 
usually faint in the Kepler bandpass (13.5 <magKep < 
16), hence the faint magnitude limit (Equation IH]) allows 


for fainter, and hence more numerous, EBs to mimic the 
transit signal than for the brighter targets in our first 
program. This higher FPR value is also because the sec¬ 
ond sample contains mainly small size planet candidates 
(super-Earth candidates) which typically have smaller 
transit depths. It is expected that the FPP increases 
with magnitude and with decreasing galactic latitude, 
and we observe such trends in figure ITS]) . Furthermore, 
the faintness of M-stars also leads to a centroid analyzes 
with lower precisions, often prevents direct-imaging, and 
consequently contributes to a higher stellar blend fre¬ 
quency. 

T he stellar parameters that we us e for this study are 
from iBatalha et ahl (1201311 . However. iHuber et al.l (|2013[1 
recently updated the stellar properties using data from 
Quarter 1-16 and different observational techniques, but 
homogeneously extract th e parameters from the Dart¬ 
mouth stellar isochrones (iDotter et al.ll2008ll . We find 
that the changes in stellar masses c orresp ond to 9% in 
average compared to iBatalha et al.l (12013 ) for the KOIs 
used for this study. There are five KOIs of our list with 
updated masses that decrease between 15% and 30% 
(KOI-104, 244, 25 2, 663, 899, and 947). As noted by 
IHuber et al.l (|2013l) . these are mainly low mass stars. Six 
stars have their updated masses that increase by more 
than 15% (KOI-13, 87, 98, 111, 446, and 1362). Since 
the mass of the primary is kept fixed in our study, we 
look at two extreme cases of mass changes (KOT252 and 
87) in order to test the effe ct of these variatio ns on our 
determination of the FPPs. IHuber et al.l (I2013D reported 
a mass for KOI-252 of 0.5 Mq, a decrease of 30%, and 
a mass for KOI-87 that is 20% larger (0.98 Mq). For 
these two objects, we compute and compare the appar¬ 
ent transit depth ratios with the old and with the new 
masses for the primary stars, and we estimate the varia¬ 
tions of these ratios for different masses of secondaries 
and tertiaries (such as presented in Fig [5]). We find 
that with a decrease in the stellar mass of 30% for a 
0.65 Mq (KOT252), the apparent transit ratios decrease 
of no more than 30% for any secondaries with masses 
lower than 1.2 Mq, and distance modulii of 5. There¬ 
fore, a primary star with a lower updated mass results 
in a lower number of low mass stars that can be rejected 
in order to satisfy the constraint from the Spitzer tran¬ 
sit depth. We compute the new FPP for KOT252.01 
considering the updated mass. We find that the FPP 
increases to 0.3%, compared to 0.21% with the previous 
mass. Inversely, increasing the mass of the primary by 
20% (0.98 Mq, KOI-87) generates an increase in the ap¬ 
parent transit ratios of less than 30% for any mass of the 
secondary that is lower than 1.2 Mq. The new FPP for 
this target goes down to 0.68%. We conclud e that the 
new stellar masses computed by IHuber et al.l (120131 ) do 
not change significantly the FPPs of the KOIs that we 
have observed with Spitzer. The FPPs presented in the 
current study are slightly underestimated for the KOIs 
fo r which the up d ated masses have decreased compared 
to IBatalha et al.l (I2013D , whereas they are slightly over¬ 
estimated for the KOIs for which the reevaluated masses 
are increased. 

Finally, the period distribution of the Spitzer sample is 
skewed towards long period candidates compared to the 
distribution obtained from the list of KOIs (see Fig. ITO)) . 
This is because the second Spitzer sample was selected 
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to focus on the super-Earth candidates that orbit in or 
close to the habitable zone of their host stars. Longer pe¬ 
riod KOIs are more difficult to characterize (less Kepler 
transit events), therefore the skewed period distribution 
would tend to overestimate the FPR extrapolated from 
the Spitzer measurements. 

8.3. Comparison with complementary studies 

The majority of the KOIs presented here have already 
been the focus of more specific studies, and have been val¬ 
idated or confirmed as true planets (see Table [3l). A sub¬ 
set of planets have constraints on their masses either from 
radial velocity measurements or fro m TTVs measured 
from the Kepler transi t lightcurves (iHolman fc Murravl 
120051: lAgol et all 120051 : iHolman et al.l l2010[) . Further- 
more, nearly half (22) of the KOIs that we selected are 
in multiple transiting systems. Eleven of these systems 
have been validated bvfR owe et al l (120141) usin g statisti¬ 
cal arguments from iLissauer et al.l (I2012ll2014ll and Ql- 
Q8 Kepler data. Most of these systems were initially 
selected from the multiples in order to understand this 
relatively unexplored class of objects. Others have com¬ 
panions that were discovered after their selection. Since 
then, we now know that candidates in multiple systems 
have a ver y high probab i lity to be planets as demon¬ 
strated by ILissauer et al.l (|20IIL 1201 2i 1201411 . However, 
for the purposes of this study, we assume that we have 
no information other than the depths and magnitudes 
measured in the Kepler and Spitzer bandpasses, and the 
direct imaging observations when available. 

In many cases the Kepler team relies on the 
BLENDER procedure to assess the planetary nature 
of candidates in a statistical manner (e.g., Keple r-9d: 
[Torres et al.ll201lL Kenler-llg: ILissauer et ^I2011L and 

Keoler-lOc: lEressin et all 120111 1. BLENDER takes into 

account the detailed morphology of the transit to reject 
as many false positive scenarios as po ssible. This ap¬ 
proach is a l so use d by other groups (e.g.. lNefs et al.l[?012l : 
iDiaz et all 1201411 . A candidate is considered validated 
when the likelihood of the signal being due to a true 
planet is much larger (by orders of magnitude) than the 
likelihood of a false positive. Many of the steps followed 
in the present analysis are inspired by the BLENDER 
approach, but are simplified and adapted to our pur¬ 
poses. In particular, we use here only the transit depth 
rather than its detailed shape to rule out blends. In 
contrast to BLENDER, we adopt a single representative 
isochrone for all stars rather than different isochrones 
for the background binary and the target based on the 
measured spectroscopic properties. This latter approxi¬ 
mation has little impact on our results. We test this for 
a sub-sample of our candidates for which we obtained 
spectroscopic reconnaissance of the main targets. The 
sub-sample comprises some candidates that were used for 
BLENDER validation and have already been published 
(see TablelSl). For this subset, we determine the mass, ra¬ 
dius, and age of the host star from a fit of the isochrones 
as described bv [Torres et al.l (j2008H . We compute a new 
FPP for each KOI in this subset. We compare this value 
to the FPP computed using the standard isochrone for 
the primary, and have checked that they are very simi¬ 
lar. This is because most of the false positive scenarios 
are mainly ruled-out from the combination of the Kepler 
photometry, the imaging and centroid information. We 


also test the robustness of our method using isochrones 
of different ages for the secondary star instead of a stan¬ 
dard isochrone. The main motivation for this test is that 
a blended unassociated triple might have a component 
that is at a very different age from the primary. We 
let the secondary isochrones range over all ages between 
0.7 to 10 Gyr, while fixing the standard isochrone for 
the primary. We find that the differences in the appar¬ 
ent transit depth ratios between Spitzer and Kepler (cf. 
Figured]) vary from 10 to 30% (in absolute) compared to 
the use of the standard isochrone for the secondary. This 
difference depends on the mass of the secondary and ter¬ 
tiary, and it is partially degenerate with distance, which 
is a free parameter here. In practice, this changes the 
blend frequency, and hence the FPP calculated in Sec¬ 
tion |3 by only a small fraction. The impact on the choice 
of isochrone ages for t he secondary has also been tested 
bv [Torres et al.l ()2011fl through the detailed study of the 
shape of transit light-curves using the BLENDER frame¬ 
work. They find that the age of the secondary does not 
change significantly their estimate of the FPP. Therefore, 
we conclude that using generic isochrones does not affect 
the overall conclusion of this paper that the FPR of our 
s ample is low. 

[Fressin et al.l ([2013D perform detailed numerical simu¬ 
lations of the Kepler targets to predict the occurrence of 
astrophysical false positives, and its dependence on spec¬ 
tral type, candidate planet size, and orbital period. They 
hnd that the global false positive rate of Kepler is 9.4%, 
peaking for giant planets (6-22 i?0) at 17.7%, reaching a 
low of 6.7% for small Neptunes (2-4i?0), and increasing 
again for Earth-size planets (0.8-1.25 R0) to 12.3%. We 
compare these findings with the sample of candidates ob¬ 
served with Spitzer that fall in a similar ove rall size range 
1.25-22R0). From the FPR estimated bv iFressin et al.l 
Ml), we can conclude that it would be extremely un¬ 
likely that we find no false positives in a random sam¬ 
ple of 51 KOIs observed with Spitzer. The difference 
betwe en our findings and the results from [Fressin et ah! 
(|2013f) is explained by the fact that our Spitzer sample 
underwent a much more stringent vetting procedure than 
typical KOIs. Indeed, the estimated FPP for our sample 
is lower than 8.8% at 3 cr level. 

We now compare our findings to the study 
of [Morton fc Johnsol ([2011h . These authors use 
the depths from the first 1235 KOIs reported 
from iBorucki et ^ (|2011h together with generic assump¬ 
tions and wi th the stellar popu lation synthesis model 
TRILEGAL (|Girardi et al.llSoOSll to derive the FPP of 
all KOIs. Their result was update d in [Morton! (120121) us - 
ing the transit depths reported bv iBatalha et al.l ( ^1311 . 
and is based on 16 months of Kepler observations. The 
main steps employed in the current paper are similar to 
their approach. One of the main differences is that we 
include color information on the transits depths, thanks 
to the Spitzer observations (depths and magnitudes). 
We also include constraints from direct imaging and 
centroids analyses when available. Another difference 
is that our work focuses on a sm a ll numb e r of K OIs, 
whereas iMorton fc Johns^ ()2011h ; iMortonI ()2012[ ) aim 
at evaluating the FPR of the complete sample of KOIs. 
Despite these differences, we find t hat our study is in 
good agreement with the study from [Morton fc JohnsonI 
(| 20 Ilf) . They predict that the FPR would be lower than 
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5% for half of the KOIs, and lower than 10% for most of 
them. Their results are consistent with the value that we 
report from our independent observational survey. Fur¬ 
thermore, we also find trends in agreement with their 
work for the FPP values as a function of Galactic lati¬ 
tude and Kepler magnitude (see Figure 116)) . 

There are other observational projects that address the 
question of the FPR in Kepler. The most comparable 
observational study to our proj ect in terms o f candi¬ 
date sample size is the work of iSanterne et all (l2012f ). 
They find a relatively high value for the FPR of 34.8% 
± 6.5% for their sample of 33 KOIs. However, there are 
several important diff erences between our ap proach and 
the one employed by ISanterne et al.l (|2012fl . The first 
difference remains in the observational method that is 
used to constrain the FPP of individual targets. Ins tead 
of using transit photometry, ISanterne et al.l (1201211 ob¬ 
tain radial velocity observations to establish the nature 
of the transiting candidates. The second important dif¬ 
ference concerns the selection of the candidates consid¬ 
ered for follow-up observations. They focus on the deep¬ 
est short-period transit signals with high SNR that Ke¬ 
pler h as detected. Unlike for our study, ISanterne et al.l 
(|2012ll restrict their targets to candidates with large tran¬ 
sit depths that are greater than 0.4%, and with short 
orbital periods lower than 25 days, and with host stars 
brighter than Kepler magnitude 14.7. This selection is 
obviously driven by instrumental capabilities. Instead, 
we select our candidates from a wider range of candidate 
sizes, orbital periods, and magnitudes that had been vet¬ 
ted by_the_CFOR Moreo ver, some of the KOIs selected 
by ISanterne et al.l (1201211 were noted as being slightly 
V-sha ped from the Kepler photometry in iBatalha et al.l 
(|2013ll . and these signals are considered as most likely 
due to EBs. We emphasize that the FPR is expected 
to be greater for larger transit depths. Therefore, one 
should expect a higher FPR for the family of giant plan¬ 
ets. 

The high FPR found in the sample 
of [Santerne_el_al] (12012|1 is consistent with the findings 
of IDemorv fc Sea^ (1201111 for which c l ose-in candidates 
were also targeted. iDemorv fc Seageil (1201111 refine the 
photometric transit light curve analysis of 115 Kepler 
giant planet candidates based on photometric data from 
quarters QO to Q2. These authors find that 14% of 
these candidates are likely false positives based on the 
detections of their secondary eclipses. 

Ground-based telescopes are also employed to examine 
the status of false-positives of a few KOIs using the same 
technique that we present here, i.e. color photometry. 
iColon et al.l (|2012ll use the GTC telescope and observe 
transits of four short-period (P<6 days) planet candi¬ 
dates. 

However, we note that the color photometric approach 
with ground-based instruments is limited to short period 
candidates for which the transit can be observed during 
the course of the night. Furthermore, the amplitude of 
the color-dependent effects for false positive detection in¬ 
creases as the two bandpasses under consideration are 
further removed in wavelength. Therefore, Spitzer is 
better adapted compared to ground-based photometric 
false-positive searches. 


We present the results from two large observational 
campaigns, which were conducted with the Spitzer Space 
Telescope, dedicated to estimating the false positive rate 
amongst a sample of Kepler candidates. We select a sub¬ 
sample of 51 candidates, spanning wide ranges in stellar, 
orbital and planetary parameter space, and we observe 
their transits with Spitzer at 4.5 /im. We measure the 
transit depths of these candidates in the Spitzer band¬ 
pass and compare them to the depths measured with 
Kepler. This technique allows us to derive the probabil¬ 
ity that a false-positive (blended eclipsing binaries) could 
mimic the transit-shape signal. We estimate that 85% of 
the KOIs from this sample have measured Kepler and 
Spitzer depths which agree at better than 3 cr level. We 
use the Spitzer observations to remove most of the red- 
dwarfs that could potentially remain as false positives. 
By combining Spitzer and follow-up observations, we es¬ 
timate that the overall false positive rate of our sample 
is estimated to be 1.3%, and lower than 8.8% at 99.7% 
of confidence. This rate implies that the vetting proce¬ 
dures of the Kepler data likely rule out a larger fraction 
of blends. Extrapolating the empirical knowledge gained 
from this small sample to the overall Kepler sample of 
candidates, we find that the overall false positive rate of 
the Kepler sample is small. In this context, at least 90% 
of the Kepler signals could be of planetary origin. 
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Figure 1. Candida te radii as function of their orbital periods (black point) for all the Kepler Object of Interests (KOIs) presented in 
IBatalha et al.l 1)201311 . Overplotted in red are the 51 KOIs that we targeted to estimate the false positive rate from Kepler and that are 
presented in this paper. We observe these 51 objects during transit in the near-infrared with Spitzer. Our selected sample spans a wide 
range of periods and sizes. 
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Figure 2. Example of a Spitzer transit light-curve observed in the IRAC band-pass at 4.5 p.m: KOI-701.03. Top panel: raw (unbinned) 
transit light-curve. The red solid line corresponds to the best fit model which includes the time and position instrumental decorrelations as 
well as the model for the planetary transit (see details in Sect. 0. Bottom panel: corrected, normalized and binned by 30 minutes transit 
light-curve with the transit best-fit plotted in red and the transit shape (with no limb-darkening) expected from the Kepler observations 
overplotted as a green line. The two models agree at 1 cr level. 































Low Kepler False-Positive Rate From Warm-Spitzer Observations 


15 





1.0005 


1.0000 


0.9995 


0.9990 L 



1.0005 

' <1 

84.01 

i : 



TT\ 

T : 

0.9995 




0.9990 





-4 

-2 0 2 

4 



1.0010 

1.0005 

1.0000 

0.9995 

0.9990 

0.9985 


1 

P 1 i ' 

103.01 I , i 



' I T I T : 


1.001 

1.000 

0.999 

0.998 


104.01 



'<^1 

<> 


-2 


-2 





Figure 3. Transit lightcurves from the Spitzer program 60028 (part-1/2). The lightcurves are obtained at 4.5 pm with the IRAC 
instrument aboard Spitzer. The data are corrected, normalized, binned in time, and combined (when multiple observations are available). 
The grey points are the measurements with their 1 a error-bars. The red solid lines correspond to the best fit model of the Spitzer data 
(unbinned) as described in Sect. [4] The transit shapes expected from the Kepler observations are overplotted as green lines. The planetary 
transit models are computed neglecting the effect from stellar limb-darkening. The name of the KOIs appear in each individual plot. 
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Figure 4. Transit lightcurves from the Spitzer program 60028 (part-2/2). Same as FigO 
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Figure 5. Transit lightcurves from the Spitzer program 80117. Same as Fig(3| 
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Figure 6. Ratios of the eclipse depth integrated in the Spitzer photometric bandpass over the eclipse depth integrated in the Kepler 
photometric bandpass as function of the mass of the tertiary star (M 3 ). These ratios are computed for an eclipsing binary stellar system 
composed of a secondary star of mass M 2 eclipsed by a tertiary stellar component (see Equation [S]). 
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Figure 7. Ratios of the dilutions integrated in the Spitzer photometric bandpass over the dilution integrated in the Kepler photometric 
bandpass as a function of the mass secondary star (M 2 ) for eclipsing binary stellar system. This system is composed of a secondary 
star eclipsed by a tertiary stellar component, that is blended with a primary star of mass M\. Two scenarios of distance modulus p are 
presented: p = 0, implying that the binary system is equidistant to the primary star (HT: hierarchical triple), and p = 5 for a background 
binary (EB: eclipsing binary) scenario. Three cases of Mi are presented for reference (see Equation. 
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Figure 8 . Same as Figure [6] but the eclipsing binary system is now blended with a primary star. These ratios are computed for an 
eclipsing binary stellar system composed of a secondary star of mass M 2 eclipsed by a tertiary stellar component that is blended with a 
primary star of solar mass (Mi=l). See Equation [4] for more details. 
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Figure 9. Color {Kepler- Spitzer) difference between the combination of the two stars (primary and secondary) as a function of the mass 
M 2 of the secondary. The calculations are presented for three mass scenarios for the primary star. The black curves represent the color 
differences for a secondary star equidistant to the primary (HT scenario) and the red curves are for a background secondary at a distance 
modulus from the primary corresponding to 1 magnitude. 
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Figu re 10. Norm a lized period distribution of the KOIs (in black) and the Spitzer sample (in red). The KOI distribution is derived 
from [Batalha et al.l <120131) . The vertical dashed lines correspond to the median value of each distribution. The two distributions are 
broadly consistent, except at the very shortest orbital periods. 
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Figure 11. Transiting planet frequency distr ibution as function of period bins (in logarithmic scale). The distribution is derived from the 
list of KOIs published in IBatalha et all 1)201311 . assuming that all the detected signals are of planetary origin. 
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Figure 12. Top: Transit depths measured from the Spitzer 4.5 pm lightcurves (red histogram) compared to their depths measured in the 
Kepler bandpass (black histogram) for 50 KOIs targeted in this program. The error bars on the Spitzer measurements correspond to 1 cr 
uncertainties. The targets are ordered by increasing transit depths measured from Kepler towards the right. Bottom: Distribution of the 
significance of the apparent transit depth differences measured between Kepler and Spitzer. The vertical dotted red line highlights the 3 a 
uncertainties; it encompasses 85% of the candidates. Two KOIs with differences greater than 6 cr do not appear on this figure (KOI-12.01 
and 13.01). The difference in the apparent transit depths is not corrected for dilution caused by the presence of a close-by companion. 
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Figure 13. Ratios of the transit depths measured from the Spitzer 4.5 pm lightcurves (-|-3 ( t) t o the transit depths measured in the 
Kepler bandpass for 50 KOIs targeted in this program. The KOIs are ordered similarly to Figure fl^ increasing transit depths from Kepler 
towards the right. 
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Figure 14. Histogram distribution of the False Positive Probability (FPP) of the Kepler candidates (KOIs) that we observed with Spitzer. 
Half of the overall sample has a FPP< 1%. The vertical dashed line shows the Scr upper limit of the FPR (8.8%) of the KOIs we present 
in this project (see Section E2J. 
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Figure 15. Top: False Positive Probability (FPP) for each Kepler candidate (KOI) that we followed-up with Spitzer as a function of the 
target star’s Galactic latitude. Bottom: FPP as a function of the Kepler magnitude. The overall FPPs increase as the Kepler magnitude 
increases and the Galactic latitude decreases. 
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Figure 16. False Positive Probability (FPP) for each Kepler candidate (KOI) that we followed-up with Spitzer as a function of the target 
star’s Galactic latitude and Kepler magnitude. The radii of the circles increase linearly as a function of the FPP value (largest circles for 
the largest FPP). The filled circle in red colors correspond to the targets for which the FPP is lower than 0.3%, which could be considered 
as validated at 3 cr level of confidence. These validated planets represent one fourth of the overall sample. The FPP clearly increases 
towards the bottom-right corner as expected. 
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Table 1 

Measurements from Spitzer observations of program ID 60028. 


KOI 

AOR 

Magnitude 

Flux 

Depthj^gp^g.^ 

Depth^p^j^er 




(mjy) 

(%) 

(%) 


12.01 

r39525120 

10.183 ±0.005 

15.179 ± 0.168 

13.01 

r39525376 

9.397 ±0.002 

31.305 ± 0.115 

42.01 

r41010688 

8.084 ±0.003 

104.939 ± 0.753 

69.01 

r41009920 

8.331 ±0.003 

83.561 ± 0.640 

69.01 

r41010432 

8.332 ± 0.004 

83.489 ± 0.768 

70.01 

r41165824 

10.855 ± 0.006 

8.176 ±0.119 

70.03 

r39437568 

10.834 ±0.012 

8.337 ±0.221 

70.03 

r41164544 

10.841 ±0.015 

8.283 ± 0.289 

72.02 

r39369984 

9.462 ± 0.007 

29.484 ± 0.449 

72.02 

r39369216 

9.476 ±0.016 

29.115 ± 1.076 

82.01 

r39420672 

9.353 ± 0.003 

32.605 ± 0.258 

82.02 

r39419904 

9.348 ± 0.003 

32.748 ±0.207 

82.02 

r39437056 

9.353 ± 0.004 

32.617 ±0.283 

84.01 

r39369472 

10.295 ±0.005 

13.689 ±0.167 

84.01 

r39370240 

10.305 ± 0.003 

13.573 ±0.100 

84.01 

r41165312 

10.312 ±0.004 

13.481 ±0.123 

94.01 

r39421440 

10.954 ±0.012 

7.463 ±0.210 

98.01 

r39421184 

10.991 ±0.004 

7.213 ±0.072 

103.01 

r39366400 

11.032 ±0.006 

6.944 ± 0.091 

103.01 

r39366144 

11.013 ±0.012 

7.067 ±0.192 

104.01 

r39420160 

10.601 ±0.004 

10.327 ±0.096 

104.01 

r41163776 

10.586 ± 0.004 

10.474 ±0.103 

137.01 

r39365888 

11.755 ±0.014 

3.567 ±0.116 

137.02 

r39369728 

11.901 ±0.023 

3.119 ±0.164 

157.03 

r41197568 

12.196 ±0.022 

2.378 ±0.117 

157.03 

r41197312 

12.194 ±0.017 

2.381 ±0.091 

244.01 

r39437312 

9.519 ±0.003 

27.997 ± 0.221 

244.02 

r39438848 

9.494 ±0.006 

28.630 ± 0.422 

244.02 

r39439104 

9.483 ± 0.002 

28.926 ±0.137 

244.02 

r41165568 

9.497 ±0.005 

28.547 ±0.310 

245.01 

r39420928 

8.084 ±0.010 

104.914 ±2.426 

245.01 

r41009664 

7.865 ± 0.002 

128.432 ± 0.526 

246.01 

r41009408 

8.538 ± 0.005 

69.088 ± 0.773 

246.01 

r41010176 

8.533 ± 0.005 

69.385 ± 0.769 

247.01 

r39368704 

11.017 ±0.010 

7.042 ±0.165 

247.01 

r39368448 

11.013 ±0.022 

7.071 ±0.350 

247.01 

r41164032 

11.034 ±0.006 

6.935 ± 0.096 

248.01 

r39370496 

12.276 ±0.128 

2.210 ±0.564 

248.01 

r41165056 

12.280 ±0.027 

2.200 ±0.132 

248.02 

r39366912 

12.298 ±0.024 

2.165 ±0.115 

248.02 

r39367168 

12.278 ±0.015 

2.204 ± 0.075 

248.02 

r39366656 

12.269 ±0.022 

2.224 ±0.109 

249.01 

r39419648 

11.016 ±0.005 

7.052 ± 0.074 

249.01 

r39421952 

11.016 ±0.005 

7.052 ± 0.074 

250.02 

r41197056 

12.521 ±0.039 

1.763 ±0.153 

250.02 

r41196800 

12.513 ±0.040 

1.776 ±0.158 

250.02 

r41196544 

12.530 ±0.028 

1.749 ±0.107 

251.01 

r39437824 

11.485 ±0.012 

4.578 ±0.127 

251.01 

r41164800 

11.461 ±0.010 

4.681 ±0.110 

252.01 

r39421696 

12.489 ± 0.029 

1.815 ±0.117 

252.01 

r41166336 

12.488 ± 0.023 

1.818 ±0.093 

253.01 

r41440256 

12.318 ±0.010 

2.125 ±0.048 

255.01 

r39420416 

11.998 ±0.015 

2.853 ± 0.094 

260.02 

r39438080 

9.320 ± 0.004 

33.626 ± 0.304 

271.01 

r39439360 

10.236 ±0.007 

14.464 ± 0.241 

271.01 

r41166080 

10.236 ± 0.009 

14.463 ± 0.293 

273.01 

r39368192 

9.953 ± 0.002 

18.763 ±0.078 

273.01 

r39367680 

9.965 ±0.012 

18.558 ± 0.496 

273.01 

r39367424 

9.959 ± 0.004 

18.661 ± 0.169 

314.01 

r44144384 

9.322 ±0.003 

33.568 ± 0.234 

365.01 

r40252928 

9.620 ± 0.004 

25.504 ±0.211 

365.01 

r40252672 

9.616 ±0.009 

25.603 ±0.528 


0.9271 

± 

0.0019 

0.688 

0.4646 

± 

0.0031 

0.228 

0.0334 

± 

0.0004 

0.018 

0.0271 

± 

0.0003 

0.007 

0.0271 

± 

0.0003 

0.009 

0.1028 

± 

0.0021 

0.057 

0.0829 

± 

0.0019 

0.055 

0.0829 

± 

0.0019 

0.062 

0.0497 

± 

0.0005 

0.032 

0.0497 

± 

0.0005 

0.033 

0.0949 

± 

0.0034 

0.056 

0.0271 

± 

0.0011 

0.038 

0.0271 

± 

0.0011 

0.018 

0.0698 

± 

0.0016 

0.110 

0.0698 

± 

0.0016 

0.016 

0.0698 

± 

0.0016 

0.070 

0.5745 

± 

0.0016 

0.478 

0.2276 

± 

0.0006 

0.170 

0.0821 

± 

0.0114 

0.069 

0.0821 

± 

0.0114 

0.075 

0.1501 

± 

0.0219 

0.081 

0.1501 

± 

0.0219 

0.154 

0.2292 

± 

0.0054 

0.168 

0.3235 

± 

0.0087 

0.371 

0.1401 

± 

0.0047 

0.046 

0.1401 

± 

0.0047 

0.138 

0.1188 

± 

0.0012 

0.125 

0.0400 

± 

0.0003 

0.041 

0.0400 

± 

0.0003 

0.011 

0.0400 

± 

0.0003 

0.043 

0.0607 

± 

0.0017 

0.051 

0.0607 

± 

0.0017 

0.049 

0.0350 

± 

0.0003 

0.056 

0.0350 

± 

0.0003 

0.035 

0.0992 

± 

0.0241 

0.056 

0.0992 

± 

0.0241 

0.105 

0.0992 

± 

0.0241 

0.089 

0.1762 

± 

0.0187 

0.078 

0.1762 

± 

0.0187 

0.135 

0.1048 

± 

0.0182 

0.086 

0.1048 

± 

0.0182 

0.087 

0.1048 

± 

0.0182 

0.218 

0.1640 

± 

0.0020 

0.174 

0.1640 

± 

0.0020 

0.106 

0.1896 

± 

0.0103 

0.108 

0.1896 

± 

0.0103 

0.280 

0.1896 

± 

0.0103 

0.105 

0.2228 

± 

0.0425 

0.280 

0.2228 

± 

0.0425 

0.315 

0.2157 

± 

0.0726 

0.152 

0.2157 

± 

0.0726 

0.187 

0.1747 

± 

0.1242 

0.094 

0.2393 

± 

0.0636 

0.143 

0.0346 

± 

0.0006 

0.039 

0.0350 

± 

0.0008 

0.013 

0.0350 

± 

0.0008 

0.005 

0.0297 

± 

0.0101 

0.029 

0.0297 

± 

0.0101 

0.041 

0.0297 

± 

0.0101 

0.029 

0.0740 

± 

0.0139 

0.101 

0.0656 

± 

0.0039 

0.086 

0.0656 

± 

0.0039 

0.078 


+0.012 

- 0.011 

+0.011 

- 0.011 

+ 0.005 

- 0.005 

+ 0.008 

- 0.006 

+ 0.008 

- 0.007 

+0.022 

- 0.023 

+ 0.016 

- 0.017 

+ 0.016 

- 0.016 

+ 0.008 

- 0.008 

+0.010 

- 0.010 

+0.010 

- 0.010 

+0.011 

- 0.011 

+0.012 

- 0.011 

+0.020 

- 0.022 

+ 0.017 

- 0.013 

+ 0.019 

- 0.018 

+ 0.017 

- 0.017 

+ 0.016 

- 0.016 

+ 0.028 

- 0.028 

+ 0.025 

- 0.026 

+ 0.039 

- 0.043 

+ 0.039 

- 0.039 

+ 0.037 

- 0.040 

+ 0.046 

- 0.049 

+ 0.040 

- 0.035 

+ 0.040 

- 0.039 

+ 0.012 

- 0.013 

+0.011 

- 0.011 

+ 0.011 

- 0.009 

+ 0.011 

- 0.012 

+ 0.009 

- 0.009 

+ 0.005 

- 0.005 

+ 0.007 

- 0.007 

+ 0.007 

- 0.007 

+ 0.038 

- 0.037 

+ 0.034 

- 0.035 

+ 0.048 

- 0.046 

+ 0.052 

- 0.051 

+ 0.049 

- 0.050 

+ 0.057 

- 0.056 

+ 0.061 

- 0.059 

+ 0.071 

- 0.077 

+ 0.032 

- 0.033 

+ 0.030 

- 0.032 

+ 0.066 

- 0.069 

+ 0.063 

- 0.067 

+ 0.068 

- 0.071 

+ 0.041 

- 0.041 

+ 0.075 

- 0.083 

+ 0.051 

- 0.052 

+ 0.047 

- 0.049 

+ 0.066 

- 0.068 

+ 0.039 

- 0.039 

+ 0.007 

- 0.007 

+ 0.015 

- 0.011 

+ 0.011 

- 0.005 

+ 0.024 

- 0.022 

+ 0.020 

- 0.022 

+ 0.019 

- 0.020 

+ 0.020 

- 0.021 

+ 0.010 

- 0.011 

+0.010 
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Table 2 

Measurements from Spitzer observations of program ID 80117. 


KOI 

AOR 

Magnitude 

Flux 

(mjy) 

Depthj(;‘ep/er 

(%) 

(%) 

87.01 

r44159488 

10.128 ±0.004 

15.965 ± 0.163 

0.0492 ± 0.0075 

0.008l°;°J2 

111.03 

r44162048 

11.206 ±0.017 

5.916 ± 0.229 

0.0616 ± 0.0012 

0.094tO;025 

174.01 

r44162560 

11.515 ±0.011 

4.453 ±0.115 

0.1039 ± 0.0026 

0.120l°;°|® 

174.01 

r44162304 

11.508 ±0.024 

4.483 ± 0.238 

0.1039 ± 0.0026 

0.098l°;°g 

446.02 

r44161024 

12.042 ± 0.023 

2.739 ±0.142 

0.0920 ± 0.0365 

0 91 O+0-060 
u.zio_Q 059 

446.02 

r44161536 

12.033 ± 0.046 

2.764 ±0.278 

0.0920 ± 0.0365 

Oil O+0-051 

U.lio_Q 054 

446.02 

r44160768 

12.047 ±0.028 

2.727 ±0.170 

0.0920 ± 0.0365 

u.uoz_0 030 

555.02 

r44162816 

12.969 ±0.136 

1.166 ±0.313 

0.0937 ± 0.0030 

f) -|qo+0.051 

U.1U9_0 053 

663.02 

r44159232 

10.751 ±0.014 

9.002 ± 0.286 

0.0693 ± 0.0083 

0.10410026 

663.02 

r44158720 

10.765 ±0.007 

8.882 ±0.142 

0.0693 ± 0.0083 

0.08310 029 

701.03 

r44163840 

11.630 ±0.016 

4.003 ±0.145 

0.0719 ±0.0108 

0.1001° 026 

711.03 

r44158976 

12.313 ±0.098 

2.135 ±0.430 

0.0698 ± 0.0029 

0.068loo|7 

817.01 

r44160512 

12.216 ±0.053 

2.334 ± 0.268 

0.1122 ±0.0576 

0 QC^Q+O-OSO 

U.UO5_0 05;L 

817.01 

r44160256 

12.211 ±0.027 

2.344 ±0.140 

0.1122 ±0.0576 

0 073+0 076 
U.U/O_0 002 

854.01 

r44164864 

12.363 ±0.027 

2.039 ±0.125 

0.1694 ±0.1658 

0 043+0 048 

U.U49_0 035 

854.01 

r44164352 

12.363 ±0.117 

2.039 ± 0.483 

0.1694 ±0.1658 

U.1DZ_0 052 

899.03 

r44165376 

11.743 ±0.011 

3.609 ± 0.092 

0.0762 ±0.0214 

U.UO/ _0.O4O 

899.03 

r44165632 

11.736 ±0.021 

3.633 ±0.171 

0.0762 ±0.0214 

0 1 37+0 044 
u.io/ _0.o45 

899.03 

r44166144 

11.748 ±0.020 

3.594 ±0.160 

0.0762 ±0.0214 

0 1 99+0.046 

u.izz_0 043 

947.01 

r44164608 

11.889 ±0.025 

3.154 ±0.175 

0.1607 ±0.0177 

nil 3+0 041 
U.iiO_0 047 

947.01 

r44165120 

11.898 ±0.021 

3.128 ±0.146 

0.1607 ±0.0177 

0 -194+0 051 
U.1Z4_0 053 

952.03 

r44159744 

12.604 ±0.037 

1.633 ±0.134 

0.1939 ±0.0543 

0 -140+0 065 

U.140_o 000 

952.03 

r44160000 

12.602 ± 0.034 

1.635 ±0.123 

0.1939 ±0.0543 

0.32410 065 

1199.01 

r44166400 

12.681 ±0.134 

1.520 ±0.403 

0.1039 ±0.0625 

0.085lO;0|4 

1361.01 

r44161280 

12.224 ±0.061 

2.316 ±0.302 

0.1419 ±0.0198 

0.1041° 0|4 
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Table 3 
FPP results. 


KOI 

O-R-S 

(*) 

Gal. long 

1 (deg) 

Gal. lat 
b (deg) 

Magnitude 

(Kepler) 

AO 

Centroid 

FPP 

(%) 

12.01 

19.7 

75.50 

7.47 

11.35 

n 

y 

0.40 

13.01 

20.7 

77.51 

16.81 

9.96 

y 

n 

0.02 

42.01 

3.1 

74.81 

18.66 

9.36 

y 

n 

0.20 

69.01 

3.4 

71.20 

10.46 

9.93 

y 

n 

0.67 

70.01 

2.0 

73.38 

14.57 

12.50 

y 

n 

0.19 

70.03 

2.0 

73.38 

14.57 

12.50 

y 

n 

0.56 

72.02 

2.8 

80.49 

18.82 

10.96 

y 

n 

0.06 

82.01 

3.7 

76.48 

20.43 

11.49 

y 

n 

0.12 

82.02 

0.2 

76.48 

20.43 

11.49 

y 

n 

0.76 

84.01 

0.9 

70.10 

10.80 

11.90 

y 

n 

2.62 

94.01 

5.7 

76.23 

7.98 

12.20 

y 

n 

0.86 

98.01 

3.6 

78.15 

16.52 

12.13 

y 

n 

0.07 

103.01 

0.4 

70.46 

9.86 

12.59 

y 

n 

2.85 

104.01 

0.8 

76.71 

20.71 

12.90 

y 

n 

0.98 

137.01 

1.5 

79.01 

8.92 

13.55 

n 

n 

25.88 

137.02 

0.9 

79.01 

8.92 

13.55 

n 

n 

22.86 

157.03 

1.7 

76.16 

8.13 

13.71 

n 

y 

6.40 

244.01 

0.5 

70.35 

14.16 

10.73 

y 

n 

0.24 

244.02 

1.3 

70.35 

14.16 

10.73 

y 

n 

0.66 

245.01 

2.4 

74.44 

17.84 

9.71 

y 

n 

0.11 

246.01 

2.1 

80.69 

15.15 

10.00 

y 

n 

0.61 

247.01 

0.5 

80.24 

19.19 

14.22 

y 

n 

0.68 

248.01 

1.7 

73.26 

10.71 

15.26 

n 

y 

1.13 

248.02 

0.3 

73.26 

10.71 

15.26 

n 

y 

1.32 

249.01 

1.1 

76.09 

17.79 

14.49 

y 

n 

2.51 

250.02 

0.5 

76.68 

17.99 

15.47 

n 

n 

7.98 

251.01 

1.1 

81.59 

10.22 

14.75 

y 

n 

3.59 

252.01 

0.6 

80.31 

15.44 

15.61 

n 

y 

0.21 

253.01 

0.6 

80.19 

18.78 

15.25 

n 

n 

6.08 

255.01 

1.3 

73.62 

14.55 

15.11 

n 

y 

1.20 

260.02 

0.6 

75.64 

14.23 

10.50 

y 

n 

1.89 

271.01 

3.0 

76.22 

17.64 

11.48 

y 

n 

0.12 

273.01 

0.2 

69.44 

13.05 

11.46 

y 

n 

0.72 

314.01 

1.1 

75.11 

13.17 

12.93 

n 

n 

9.54 

365.01 

1.9 

83.14 

11.64 

11.19 

y 

n 

0.92 


87.01 

2.9 

79.09 

15.79 

11.66 

y 

n 

0.70 

111.03 

1.2 

73.19 

14.57 

12.60 

n 

n 

5.87 

174.01 

0.1 

81.57 

11.31 

13.78 

n 

n 

14.09 

446.02 

0.4 

69.82 

16.13 

14.43 

n 

n 

8.74 

555.02 

0.2 

73.86 

10.31 

14.76 

n 

n 

34.51 

663.02 

1.2 

72.20 

16.03 

13.51 

y 

n 

6.19 

701.03 

1.0 

75.04 

18.69 

13.73 

n 

n 

12.14 

711.03 

0.1 

79.45 

11.29 

13.97 

n 

n 

22.19 

817.01 

0.6 

69.88 

16.29 

15.41 

n 

n 

10.44 

854.01 

0.4 

73.45 

13.13 

15.85 

n 

n 

21.25 

899.03 

0.7 

77.66 

9.07 

15.23 

n 

n 

16.97 

947.01 

1.1 

78.48 

13.56 

15.19 

n 

n 

12.40 

952.03 

0.6 

80.54 

9.96 

15.80 

n 

n 

42.18 

1199.01 

0.2 

72.28 

8.97 

14.89 

n 

n 

19.32 

1361.01 

0.7 

76.01 

9.58 

14.99 

y 

n 

16.97 


comments 


Kepler-13b Bani^^^aL 
Kepler-410A b 
Kepler-93b 
Kepler-20b 

Kepler-20c Gautie^e^aL 
Kepler- 10c FVessh^^aL 
Kepler-102e ^arc^^^al^ 
Kepler-102d ^arc^^^al^ 
Kepler-19b Ballar^e^al. 
Kepler-89d Weis^^aL 
KeDler-14b Buchhave et al. 

TTVs - 



Kepler-94b 
Kepler-18c 

Kepler- 18 d^^^^^^^^^r 

Kepler-1 lb ^issaue^^ay ^OH) 
Kepler-25c Steffei^^aL 
Kepler-25b^^^^^^^^T 
Kepler-37b ^arcl^^^aL ^^^3) 
Kepler-68b ^illilan^^^ay 


Kepler-49b Steffer^^aL ^0^) 
Kepler-49c^^^^^^^^y 


Kepler-26c^teffer^^aL ^0^) 
Kepler-125b ^owe_et_^ ^201^ 


Kepler- 126d^ow^^L (^0^ 
Kepler-127d^owe_et_^ (201^ 

Kepler-138c^£w^e^aL 


Kepler-22b£oruck^^aL ^0^) 
Kepler-104d^ow^^^aL 

Kepler-157b^ow^^^aL (2^^) 

Kepler-205c^owe_et_^ (201^ 
Kepler-62e ^oruck^^aL 

Kepler-236c^ow^e^aL (2^^^ 

Kepler-249d^o^_et_^ (201^ 

Kepler-32b ^]abr^ck^^^al^ {2Q\2^ 

Kepler-61b^allar^^^^L 


Difference between the Spitzer and Kepler apparent transit depths in unit of cr from combined measurements. The 
differences are not corrected for dilution caused by the presence of a close-by stellar companions. 

These columns indicate whether Information from the Adaptive-Optic imaging (AO) followup and Kepler centroid 
analysis are available: ”y” means that information is available and used in the study, ”no” means that no information on 
AO and centroid have been used. 
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